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I = 2 7T7T scattering using G-parity boundary condition * 
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To make the irn state with non-zero relative momentum as the leading exponential, we impose anti-periodic 
boundary condition on the pion, which is implemented by imposing G-parity or H-parity on the quark fields at the 
boundary. With this, we calculate the 1 = 2 nn phase shift from lattice simulation by using Liischer's formula. 
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1. Introduction 

Lattice gauge theory provides a way to inves- 
tigate low energy physics of QCD, which cannot 
be done using any perturbative method. One of 
the interesting physical quantities is related to 
7T7T system. The K — > tttt WME which violates 
CP symmetry is of particular interest. Since lat- 
tice calculations easily extract ground state, the 
need to generate the final 7T7t state with non- 
trivial relative momentum is a serious difficulty. 
We proposed to use G-parity boundary conditions 
to overcome this difficulty [2]. In the present pa- 
per, we have implemented this G-parity idea and 
a new H-parity boundary condition in numerical 
simulations and calculated the 1 = 2, tttt phase 
shift. 

2. G parity boundary condition 

Since the G-parity operation on a pion gives 

G\i^ >= -fa* > , G\ir° >= — |tt >, (1) 

by applying this operation on the boundary, we 
can impose anti-periodic boundary condition on 
pion. To implement this condition on lattice, we 
have to use the G-parity operation on the quark 
fields: 



(2) 



In actual calculation, we impose this boundary 
condition only in the z-direction so that we have a 
pion with non-zero z-momentum. Because at the 
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boundary there are terms such as ipip and t^ip, it 
requires some special care to implement [3]. First, 
we have to impose a charge-conjugate bound- 
ary condition on the gauge field to keep gauge 
invariance and we have to virtually double the 
box size for the Dirac operator inversion. Since 
isospin has an important role in the two pion sys- 
tem, it is worth noting that G-parity commutes 
with isospin, which means that under this un- 
usual boundary condition isospin is still a good 
quantum number. 

3. H parity boundary condition 

An easier way to impose anti-periodic bound- 
ary conditions on a pion is to apply the following 
operation on the quark fields, 



(3) 




We will call this operation H-parity 
operation on the pions will be, 



H\tt u >- 



k° > 



Then, the 



(4) 



Under this boundary condition, isospin is not a 
good quantum number anymore, but I z is still 
good. Since we know that the I z = 2, 7T7t state is 
composed of two 7r + , this state must have non- 
zero relative momentum. This is not true for the 
7 = state. So the utility of this boundary con- 
dition is more limited than the G-parity bound- 
ary condition. However it has the advantage that 
it doesn't require any modification of the gauge 
field boundary condition. This allows us to use 
existing lattices, including dynamical ones. 
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Figure 1. 1/a = 0.978^!^GeV, /3 =5.7, 8 3 x 
32, L s = 10, Wilson Gauge Action, Domain Wall 
Fermions (DWF), H-parity boundary condition. 
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Figure 2. 1/a = 0.978(l T 4^GeV, (3 =5.7, 8 3 x 32, 
L S =10, Wilson Gauge Action, DWF, H-parity. 

4. Results 

4.1. Single pion 

We first investigate the properties of the one- 
pion system. We expect to find a one-pion state 
with momentum unlike the conventional =r. 
Figure 2] shows a graph of energy versus pion 
mass. As expected, the single pion state with 
smallest energy has E{m v ) = ^Jm\ + sin 2 (j). 

4.2. Two pions 

Figure |2 shows the effective mass plot for the 
two pion state. We have a very nice plateau in the 
time range from 3 to 11, but after that we have 
suspicious fall-off. This fall-off can be explained 
by considering Fig. [3] A pion created at t = can 
propagate in either direction. Because we have 
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Figure 3. One pion "state exists all time. 

two particles, there is a state in which the two 
particles propagate in opposite directions. The 
correlation function for this state is 

G(t) « e -EAT-*) e -E w t = £ -E„T = e -2E^ (g) 

here E„ means the energy of one pion with mo- 
mentum. This is just a constant, and will be dom- 
inant near t = ^ because the energy of the 1 = 2, 
two-pion state is bigger than 2E V and can cause 
the abrupt fall-off in this region. We confirmed 
this idea by fitting the effective mass plot with the 
function "cosh + const? , the solid line in Fig. 

Figure 0] shows the same two-pion state effec- 
tive mass plot for the G-parity boundary condi- 
tion. This G-parity effective mass plot is quite 
different from the one for H-parity. Instead of 
an abrupt fall off, it has gradual decrease. It 
even looks like it has two plateaus. A simula- 
tion with more time slices (N t — 48) also shown 
in Fig. 2] demonstrates this two-plateau struc- 
ture. Since the spatial lattice volume was small 
(~ (1.7/m) 3 ) for this calculation, we guessed 
that it might be a finite volume effect and per- 
formed the same simulation with a bigger volume. 
Figure [5] shows the result for a spatial volume 
« 1.7/m x 1.7/m x 3.4/m. 

We notice that the gradual decrease of G-parity 
plot has disappeared and it is almost identical to 
that seen with H-parity. Therefore, we can con- 
clude that it is a finite volume effect that caused 
the two plateau behaviour. This might mean that 
we have discovered a new finite volume qqqq state. 
After becoming convinced from these tests that 
we have a two-pion state with non-zero relative 
momentum, we extended this simulation and cal- 
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Figure 4. 1/a = 0.978(l^eV, (3 
L s =10, Wilson Gauge Action, DWF 
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Figure 5. 1/a = 0.978 (lljGev, (3 =5.7, 8 x 8 x 
16 x 32, L S =10, Wilson Gauge Action, DWF 



culated the I — 2 tttt phase shift. Since we are 
using Luschcr's formalism this is nothing but 
spectroscopy. Figure HO shows our phase shift cal- 
culation including CP-PACS [J and experimental 
results. The following are our 5^ simulation pa- 
rameters (1/a is in GeV and iV t =32): 



Domain Wall Fermions, L 
p w 250MeV 

G 8 2 x 16 

H 8 2 x 16 
p w 450MeV 

H 8 3 

H 16 3 



=10, M 5 =1.65 



1/a 

1/a : 

1/a 

1/a : 



0.978(14) 
0.978(14) 

0.978(14) 
1.98(3) 



Wilson 91conf 
Wilson 172conf 

Wilson 269conf 
DBW2 156conf 
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5. Conclusion 

We have tested the idea of imposing anti- 
periodic boundary conditions on the pion by ap- 
plying a G-parity or H-parity operation on the 
quark field at the boundary. For the one pion 
case, we found that the energy of the particle is 
given by \Jm\ + sin 2 (j) as expected. We can 
achieve a two-pion state with non-zero momen- 
tum, and from this relative momentum, we can 
calculate the 1 = 2, tttt phase shift. Since the 
H-parity boundary condition can be applied to 
existing lattices, it will be more convenient than 
G-parity. In particular, we plan to use the H- 
parity boundary condition for a AI = | K — > tttt 
calculation on existing lattices. However, for the 
1 = 0, tttt state only the G-parity boundary con- 
dition with dynamical lattices will work. Since 
the G-parity boundary condition is more vulner- 
able to finite volume effects, e.g. the qqqq state 
which we need to study further, it may require 
more resources. 
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